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We present a study of triple Higgs boson (3H) production at the International Linear Collider 
(ILC) within the general Two-Higgs-Doublet Model (2HDM). We compute the production cross- 
sections at the leading-order for the 3H final states and find values up to <7 0.1 pb. This result 
represents a large enhancement with respect to the corresponding MSSM cross-sections, which 
stay typically at the level of a ~ 10^^ pb or less. Furthermore, since the 3H cross-sections in 
the general 2HDM can be of the order of the double Higgs production cross-sections, such 3H 
processes could be a competitive (if not the dominant) mechanism for Higgs boson production 
at the ILC. In practice, these 3H events could be identified through the tagging of 6 heavy-quark 
jet final states and, in this case, they would provide strong evidence of an extended Higgs boson 
sector - likely of non-supersymmetric nature. 
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1. Introduction 

Double Higgs boson (2H) production in a linear collider e+e- h°A°;H^^A°;H+H-, has 
been investigated in great detail in the literature, although mainly in the Minimal Supersymmetric 
Standard Model (MSSM)[|l]]. However, a tree-level analysis of these pairwise-produced uncon- 
ventional Higgs bosons is most likely insufficient to disclose their true nature. It is for this reason 
that a dedicated work on radiative correction calculations on Higgs boson production has been 
undertaken in various models |Q, Besides, a second interesting mechanism for studying the 
Higgs boson properties is the triple Higgs (3H) production, which carries essential information 
to trace back the ultimate structure of the corresponding Higgs potential. This kind of processes 
have been profusely studied also within the MSSM, although in this case the cross-sections are 
rather meager [Q]. Our main purpose here is to study the trilinear coupling HHH in the general 
Two-Higgs-Doublet Model (2HDM) by focusing on the 3H final states produced at the ILC |j5|] : 

e^e' ^H+H-hi, e^e'^hihiA^, e^e'^h^H^A^, {hi = h^ ,A^). (1.1) 

Let us recall that the general 2HDM [|6|] is obtained by canonically extending the SM Higgs sec- 
tor with a second SU{2)l doublet with weak hypercharge F = 1, so that it contains 4 complex 
scalar fields. The free parameters A,- in the general, CP-conserving, 2HDM potential can be finally 
expressed in terms of the masses of the physical Higgs particles (M^o,MffO,M^o,MH±), tanjS (the 
ratio of the two VEV's (H^) giving masses to the up- and down-like quarks) and the mixing angle 
a between the two CP-even states. There remains, however, the coupling A5, which cannot be 
absorbed in the previous quantities. Following [^, we set A5 = = 2\/2Gf M^o- This condition 
allows to keep closer to the MSSM structure of the Higgs sector. Therefore we end up with 6 free 
parameters, to wit: (M/,o,M^o,M40,M//±,tana,tanj8) . Furthermore, to ensure the absence of tree- 
level flavor changing neutral currents (FCNC), two main 2HDM scenarios arise: 1) type I 2HDM, 
in which one Higgs doublet couples only to down-like quarks, whereas the other doublet does not 
couple to any quark; 2) type II 2HDM, where one doublet couples only to down-like quarks and 
the other doublet to up-like quarks. The MSSM Higgs sector is actually a type II one, but of a very 
restricted sort (enforced by SUSY invariance) 

Further constraints must be imposed to attest that the SM behavior is sufficiently well repro- 
duced up to the energies explored so far, namely: /) the perturbativity and unitarity bounds; //) the 
approximate SU{2) custodial symmetry, i.e. \5p2hdm\ < 10^^ [^, and ///) consistency with the 
low-energy radiative B-meson decay (which entails Mh± > 350 GeV for tanjS > 1 [^, unless we 
consider type-I 2HDM). We refer the reader to ^ for further details, in particular for the full list 
of trilinear couplings within the general 2HDM that are relevant for the present calculation. 

2. Double and triple Higgs boson production in the 2HDM 

In what follows we shall mainly discuss the leading-order 3H production at the ILC within 
the general 2HDM and compare briefly with the 2H production processes. Throughout the present 
work, we have made used of the standard computational packages [p^]. 

We begin by reporting on the tree-level results for the 2H production within the 2HDM, specif- 
icaUy on the CP-conserving channels e^e^ — > h^A^;H^A^;H^H^. Let us firstly define two diffe- 
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rent scenarios: (a) light, and (b) heavy, Higgs boson masses (see Table [T]). We wish to use set I and 
III for the study of 2H production, while set II and set III for the study of 3H production. 





Set I 


Set II 


Set III 


M,,o {GeV) 


100 


100 


200 


Mh± (GeV) 


120 


120 


350 


Mfjo (GeV) 


150 


150 


250 


M^o (GeV) 


140 


300 


340 



Table 1: Sets I, II and III of light and heavy Higgs boson masses in the 2HDM. Sets I and III are used for 
2H production and Sets II and III for 3H production. 

We have found [Q] that, in the light Higgs boson mass regime, the 2H production rates within 
the 2HDM are substantial (a few thousands events per 100 fb~' of integrated luminosity), with 
maximum values reaching a{e^ — > H ) ~ 0. 1 pb, whereas for heavier Higgs masses (Set 
III) the optimal production rates lie around one order of magnitude below. Nevertheless, even 
in these less favored scenarios, the predicted rates are still quite sizeable within the clean ILC 
environment. Worth noting is that the predicted cross-sections for the same 2H processes within 
the MSSM yield a ~ 10^^ pb and are thus comparable to the 2HDM values for similar masses. 
This is a reflex of the fact that formally the Higgs-Higgs-gauge boson couplings do not differ from 
one model to the other. Therefore, we conclude that sizeable rates of non-standard 2H production 
can be achieved at the ILC for both SUSY and non-SUSY extended Higgs sectors, which implies 
that both scenarios are hard to distinguish at the tree-level. A clear separation of them can only 
be accomplished through the detailed study of radiative corrections to 2H production in both the 
MSSM [§] and the 2HDM [|]. 

Let us now discuss the case of the triple Higgs boson production (cf. Eq. ( |l.l[ )) within the gen- 
eral 2HDM (see Ref. [§] for further details and discussions). The basic result here is very different 
from the 2H case sketched above, in the sense that the 3H cross-sections for the general 2HDM 
may carry large enhancements at the tree-level which are not present in the MSSM. It means that 
the 3H channels in the general 2HDM could be, in contrast to the 2H ones, truly distinctive already 
at the leading-order. The dynamical reason for this stems from the structure of the trilinear Higgs 
boson couplings HHH in the general 2HDM. In contrast to the Yukawa coupling with fermions, 
these trilinear couplings do not depend on whether we are in type I or type II models, and can be 
largely enhanced in certain regions of the parameter space (see Table 1 of [Q]). Such enhancements 
are not possible in the MSSM case, owing to the purely gauge nature of these couplings in the 
SUSY case. 

The numerical analysis fully corroborates our expectations. In Fig. |l| we have plotted the 3H 
production cross-sections within the 2HDM. Again, two different mass regimes (sets II and III in 
Table |l|) are considered. This is so because of the constraints imposed by the radiative B-meson 
decays on the charged Higgs boson masses for type II models (M//± > 350 GeV) |^. Lighter 
Higgs boson regimes (allowed for type I models) entail optimal cross-sections at the level of ~ 0.1 
pb or more, therefore implying promising rates of at least 10^ events per 100 fb^' of integrated 
luminosity. The corresponding results for the heavy Higgs boson scenario lie around 1 — 2 orders 
of magnitude below those for the light Higgs boson scenario, and they are attained at higher values 
of ^/s. The results for set III (suited for type II models) translate into rates of 0(10^ — 10^) events 
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per 100 fb^' of integrated luminosity, which may still allow some comfortable detection of the 
signal in the clean ILC environment. The remaining 3H channels (among the 7 CP-conserving 
triple final states in (|1.1[)) provide smaller production rates in optimal conditions 
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Figure 1: Total cross-section CT(pb) for the triple Higgs boson production processes e+e^ iP , e+e^ 

//+ H^, e+ li^ ifA^ and e+ /f" A° in the general 2HDM as a function of -/s and for different values 
of tan /3 . In each case the label of the process and the choice (Set II or Set III) of Higgs boson masses used for the 
calculation is indicated, see Table [l| 

In order to compare the 2HDM results with the corresponding MSSM values, we have system- 
atically searched for the (allowed) regions across the MSSM parameter space where the optimal 
values for the cross-sections are attained. Besides the exception triggered by the resonant process 
e+e^ — > hPh^A^, whose maximum cross-section reads cj(v^ = 1 TeV) ~ 10^^ pb, the remaining 
(7(3//) are very small, namely they may reach o{^/s = 1 TeV) ^ 10^^ pb at most. We can thus 
assert that most of the 3H cross-sections in the MSSM are really tiny and, hence, very difficult to 
detect in practice. Let us finally stress that the extremely clean environment of the ILC should allow 
a relatively comfortable tagging of the three Higgs boson events for the typical 3H cross-sections 
that we have obtained in the general 2HDM case. 



3. Discussion and conclusions 

We have computed the leading-order cross-sections for triple Higgs boson production in a 
e^e^ linear collider - cf. Eq. ( |l.l| ) - within the general 2HDM. We have shown that the production 
cross-sections may comfortably reach O.lpb within type I 2HDM. Such optimal rates are achieved 
for sufficiently large (or small) values of tanjS (tanjS > 20,tanj8 < 0.1). Moreover, in certain 
regions of the parameter space, and for the most favorable processes (such as —>■ H h^) 
the cross-sections can be pushed up to ~ Ipb. In spite of the fact that 3H production does not 
involve any kind of Higgs-fermion interaction, and hence the predicted rates should not depend 
on whether type I or type II 2HDM is considered, low-energy B-meson physics puts stringent 
constraints on M//± in the type II case. Such scenarios with relatively heavy Higgs boson masses 
render maximum cross-sections that are roughly 10 times smaller, i.e. of order of 0.01 pb. The 
smaller number of events falls nonetheless in the range of 10^ per 100 fb^^ of integrated luminosity. 
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Remarkably, for both type I and type n 2HDM models the maximum cross-sections lie far above 
their MSSM counterparts (which typically remain at the tiny level of a ~ 10~^pb). 

Worth stressing is also the fact that, for the general 2HDM, the maximum 3H production rates 
are fully comparable, and even larger (for type I models), than the optimal values achieved by the 
2H processes (a ~ 0.1 pb). Such a result can be traced back to the potential enhancement of the 
Higgs boson self-interactions (HHH) within the 2HDM, which cannot be realized in the MSSM 
case because of the SUSY invariance of the interactions. We conclude that the 3H final states in the 
general 2HDM can be a competitive, if not the dominant, Higgs boson production mechanism at 
the ILC. Owing to the extremely clean ILC environment, we expect that the leading 3H signatures 
could hardly be missed. The latter could be quite spectacular since they should reveal in the form 
of 6 heavy-quark jet final states. If a few, well identified, events of this kind would be detected, it 
should hint strongly at (non-SUSY) Higgs boson physics beyond the SM 
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